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We analyze the unpolarized and polarized electron-induced proton- knockout reactions on O in 
different kinematical settings using two theoretical approaches. The first one is based on a relativistic 
mean-field distorted-wave description of the bound and scattering states of the proton, including a 
fully relativistic electromagnetic current operator. The second approach adopts the same current 
operator, but describes the proton properties consistently on the basis of microscopic calculations of 
the self-energy in 16 O below the Fermi energy and final-state damping in nuclear matter above the 
Fermi energy, using the same realistic short-range and tensor correlations. Good agreement with 
all unpolarized and polarized data is obtained at low and high Q 2 by using the same spectroscopic 
factors fixed by the low-Q 2 analysis, indicating that a high degree of internal consistency has been 
reached. 
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A long series of high-precision experiments on sev- 
eral nuclei [Q, |2| ||, ^, 5j have generated a well es- 
tablished tradition which singles out exclusive (e, e'p) 
knockout reactions as the primary tool to explore the 
single-particle aspects of the nucleus. The experimen- 
tal analysis has focused on the missing energy spec- 
trum of the nuclear response, assigning specific quan- 
tum numbers and spectroscopic factors to the various 
peaks corresponding to orbitals close to the Fermi en- 
ergy. In addition, the missing momentum dependence 
of these spectra has been studied, stimulating for ex- 
ample the exploration of the high-momentum compo- 
nents induced by nucleon-nucleon correlations inside 
nuclei || 0] . The theoretical description of these reac- 
tions have usually been performed in the framework 
of the nonrelativistic distorted- wave impulse approxi- 
mation (DWIA), including the Coulomb distortion of 
the electron and proton waves due the presence of the 
nuclear field g, [|, [u^. This approach was able 
to describe to a high degree of accuracy the shape of 
the experimental momentum distribution for several 
nuclei in a wide range of different kinematics ]9|, 10 . 
However, a systematic rescaling of the normalization 
of the bound state, interpreted as the spectroscopic 
factor for the corresponding level, had to be applied 
in order to reproduce the magnitude of the experi- 
mental distribution ^ [ll|. This systematic deviation 
from the mean-field expectations has clearly identified 
the limits of this approximation. In fact, nowadays a 
clear picture has emerged in which a considerable mix- 
ing between single-hole states and more complicated 
configurations results in a fragmentation of the single- 
particle strength in several peaks around and beyond 



the Fermi surface. A further depletion of the single- 
particle strength is induced by short-range and ten- 
sor correlations between nucleon pairs in the ground 
state 

More recent (e, e'p) experiments have been carried 
out at the Jefferson Laboratory (JLAB) |l4|, [l5) at 
higher momentum transfer Q 2 and with increased 
statistics such that the fully differential cross section is 
now directly available. The new kinematic domain re- 
quired a substantial upgrade of several theoretical in- 
gredients in order to incorporate all possible relativis- 
tic effects. Models based on relativistic DWIA (RD- 
WIA) have been developed, where the Dirac equation 
is solved directly for the nucleon bound and scattering 
states |ll| 17, |ll| |l9| or, equivalently, a Schrodinger- 
like equation is solved and the spinor distortion by 
the Dirac scalar and vector potentials is incorporated 
in an effective current operator in the so-called effec- 
tive Pauli reduction |2Cj . A successful description 
of the data has been achieved, but slightly different 
spectroscopic factors are deduced, because the rela- 
tivistic optical potentials in general give a stronger 
residual final-state interaction (TSI) than the corre- 
sponding nonrelativistic ones [[18| pifl . Moreover, the 
limits of validity of the older DWIA analysis versus 
RDWIA were not always properly explored, as dis- 
cussed in Ref. p2| , resulting, for example, in a certain 
degree of ambiguity for the spectroscopic factors ex- 
tracted at low energy. 

Despite several sources of theoretical uncertainties 
(different equivalent potentials for FSI, relativistic ef- 
fects on both FSI and spectroscopic factors, off-shell 
effects...), a consistent microscopic treatment of the 
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(e, e'p) reaction mechanism at different kinematics is 
highly desirable. Results for a first attempt towards 
this goal were recently obtained in Ref. [£3| (see also 
Ref. concerning the treatment of FS1), where 
a successful analysis of low- and high-Q 2 data was 
performed using identical spectroscopic factors which 
were deduced at low Q 2 . In the present paper, this 
analysis is extended to the corresponding JLAB ex- 
periment with polarization, as reported in Ref. |l5| . 
The results will then be compared with those obtained 
in the RDWIA approach of Ref. [^2|, where a con- 
sistent description of low- and high-energy data was 
generated and a careful analysis of the limits of the 
nonrelativistic DWIA was carried out. The sensitivity 
to different off-shell prescriptions for the electromag- 
netic current operator will be also discussed |2{| , but 
the difference between spectroscopic factors obtained 
by nonrelativistic and relativistic analyses remains un- 
solved and its discussion is beyond the scope of this 
paper. 

The basic ingredient of the calculation is the tran- 
sition amplitude (omitting spin degrees of freedom for 
simplicity) || 

J^,q,f N ,E R ) = J dp dp' X^XjP') 

x J^(p,p',q,iv) 4> ERn (p) [Z n (E R p , (1) 

where q, to are the momentum and energy transferred 
to the target (Q 2 = q 2 — u> 2 ) and p' N is the knocked- 
out nucleon momentum, leaving the residual nucleus 
in a well-defined state with energy Er and quantum 
numbers n. The function </>£ B „ describes the over- 
lap between the exact A— body initial state and the 
residual (A — 1)— body state induced by producing a 



hole; x 



(-) 



describes the same kind of overlap when 



producing the hole in the exact A-body final state ||] . 
The norm of <f> E n is 1 and Z n (En) is the spectro- 
scopic factor associated with the removal process, i.e., 
it corresponds to the probability that the residual nu- 
cleus can indeed be considered as the target nucleus 
with a hole. The boundary conditions of the eigen- 
value problem for v^, L are those of an incoming 
wave. 

In the RDWIA of Refs. ||, ||, £~ „ is replaced by 
the solution of a Dirac equation [E6f deduced in the 
context of a relativistic mean-field theory that satis- 
factorily reproduces global and single-particle proper- 
ties of several nuclei |27| . For the scattering states the 
effective Pauli reduction is applied. The Darwin non- 
locality factor, that contains the effect of the negative- 
energy components of the spinor, is reabsorbed in the 
current operator, which becomes an effective relativis- 
tic one-body operator depending on the Dirac scalar 
and vector potentials 19, 2(J, as well as on the chosen 
off-shell prescription (ecl, cc2, or cc3) |^5|, ^8| . The 
function ~ yL ' becomes a two-component 

spinor which solves the Schrodingcr equation with 



the equivalent central and spin-orbit potentials ex- 
pressed in terms of the original Dirac scalar and vector 
ones |29f . 

In Ref. [^3| the transition amplitude is evaluated 
by systematically applying the effective Pauli reduc- 
tion to both the initial and final Dirac spinors, de- 
termining the relevant integrals in momentum space 
thus avoiding any effective momentum approximation 
(EMA) |2C[| . The current operator displays the same 
features as in the RDWIA discussed above, i.e., it is an 
effective one-body relativistic operator depending on 
the Dirac scalar and vector potentials. In Ref. |23| 
only the ecl off-shell prescription has been consid- 
ered for compatibility with older low-energy data anal- 
yses Q. The scattering state of the (very ener- 
getic) proton is described in the eikonal approxima- 
tion by a uniformly damped plane wave, or, equiv- 
alently, by a plane wave with a complex momentum 
Pf — PN+i-Pi H3' EH' imaginary part pi is micro- 
scopically justified by linking the proton absorption to 
the same process taking place in nuclear matter and by 
calculating the nucleon self-energy in a self-consistent 
manner with realistic short-range and tensor correla- 
tions p3| , |33| . The observed damping is also in agree- 
ment with experimental expectations in different kine- 
matic domains p4j , however, embedding the proton in 
nuclear matter prevents the inclusion of spin-orbit ef- 
fects; therefore, the corresponding Darwin nonlocality 
factor for the final state is just 1. The function (j) E n 
is obtained from p-shell quasihole states deduced from 
the nucleon self-energy calculated for 16 O using real- 
istic short-range and tensor correlations [Q. 

The hadronic tensor of the reaction, W^" , involves 
an average over initial states and a sum over the un- 
detected final states of bilinear products of the scat- 
tering amplitude (Q). The differential cross section for 
the (e, e'p) reaction, with initial beam helicity h and 
proton polarization component s, becomes Q 

day, s e 4 1 

If- = L W^ v 

dp'Jp' 16^ 2 Q±p eP ' e ^ 
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167r 2 Q 4 p e p[ 

da° 1 



dp£dp' N 2 



i J f 



l + P-s + h(A + P' -s) 



(2) 



where p e ,p' e are the initial and final electron momenta 
and L^y is the lepton tensor. The coefficients of 
the linear expansion are the induced polarization P. 
the electron analyzing power A, and the polarization 
transfer coefficient P' . The reference frame in the 
polarimeter is formed by the direction of p' N (L com- 
ponent), the direction of q x p' N (N component) and 
N x L (T component). In coplanar kinematics, as is 
the case for the E89033 experiment at JLAB |l5| , only 
pN p/L anc j p/T surv j va When summing over the 
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recoil proton polarization and the beam helicity, the 
usual unpolarized cross section da° is recovered. 
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FIG. 1: Upper panel: cross section for the le O(e, e'p) 15 N 
reaction at E p = 90 MeV constant proton e nerg y in the 
center-of-mass system in parallel kinematics [B5|. Lower 
panel: cross section for the same reaction but at Q 2 = 0.8 
(GeV/c) 2 in perpendicular kinematics p^ |. Data for the 
pi state have been multiplied by 40 and 20, respectively. 
Solid lines show the results when using the quasi-hole spec- 
tral function for the bound state (see text) with spectro- 
scopic factors Z pl /2 = 0.644 and Z pS /2 = 0.537 in both 
panels p3| |3c| ] . Dashed, dot-dashed, and dotted lines rep- 
resent the result of the RDWIA approach with ccl, cc2, 
cc3 off-shell prescriptions, respectively (see text). All the 
RDWIA curves in both panels have been rescaled by the 



spectroscopic factors Z vl / 2 = 0.708 and Z t 
obtained by a y 2 
current. 



J p3/2 



0.602, 



fit to the data of Ref. |35[| using the cc3 



In Fig. [I] we first reconsider the unpolarized 
16 0(e,e'p) reaction leading to the ground state and 
the first excited state of 15 N with p^ and p| quan- 
tum numbers, respectively. In the upper panel, data 
have been collected in parallel kinematics {p' N \\ q) at 
a constant proton energy of 90 MeV in the center-of- 
mass system They are presented in the form of 
the reduced cross section 



n(p m ,E m ) = 



1 



da° 



dp' e dp' N Ka 



(3) 



as a function of the missing momentum p m = p' N — q 
at the considered missing energy E m , where K is a 
suitable kinematic factor and a ep is the elementary 
(half off-shell) electron- proton cross section p8|. For 



ease of viewing, the results for the transition to the p\ 
ground state have been multiplied by 40. The solid 
lines refer to the calculations employing the p-shell 
quasihole states for 16 O in a nonrelativistic frame- 
work, as discussed in Ref. pfl. The spectroscopic fac- 



tors extracted from the data are Z, 



pl/2 



0.644 and 



Zp3/2 = 0.537, respectively. The dashed lines show 
the results of the RDWIA analysis with the same ccl 
off-shell prescription; dot-dashed and dotted lines in- 
dicate the results when using the cc2 and cc3 recipes, 
respectively. Hence, the comparison among dashed, 
dot-dashed, and dotted lines shows the evolution in p m 
of the theoretical uncertainty related to offshcllness at 
this kinematics. Relativistic mean-field bound states 
are obtained by solving Hartree-Bogoliubov equations 
with finite-range interactions [p6[ . The proton scat- 
tering wave is deduced from relativistically equiva- 
lent energy-dependent optical potentials p9[ . The 
resulting spectroscopic factors, Z p u 2 — 0.708 and 
Zp3/2 — 0.602, have been obtained by a x 2 fit us- 
ing the cc3 current, which gives an overall better de- 
scription of the (e, e'p) observables, particularly for 
the left-right asymmetry. In the lower panel, the same 
reaction is considered at constant (q, u>) with Q 2 = 0.8 
(GeV/c) 2 @. The data now refer to the fully differ- 
ential unpolarized cross section da" , avoiding any am- 
biguity in modeling the off-shell behavior of <r ep J28| . 
Again the solid lines refer to the calculation employing 
the p-shell quasihole states but with an effective rel- 
ativistic current operator and an eikonal microscopic 
description of FSI as discussed above (see Ref. (23) for 
further details). The dashed, dot-dashed, and dot- 
ted lines still refer to the RDWIA analysis with the 
ccl, cc2, and cc3 off-shell prescriptions for the elec- 
tromagnetic current, respectively. The p^ results are 
multiplied by a factor 20. The theoretical curves are 
rescaled by the same spectroscopic factors as in the 
upper panel. The agreement with the data remains 
very good also in this case. This confirms the in- 
ternal consistency of the two approaches, since the 
spectroscopic factors correspond to a nuclear prop- 
erty that must be independent of the probe scale Q 2 . 
Incidentally, we remark that an extraction of spec- 
troscopic factors directly from the data of Ref. JTJ] 
most likely produces ambiguous results, due to the 
small number of data points (8 only) of the experi- 
ment. We tried such an extraction using the RDWIA 
curves, but the fits had very high \ 2 P er degree of 
freedom and gave p\ rescaling coefficients which are 
systematically bigger than the p^ ones, contrary to 
any reasonable expectation (see also Ref. [Q). The 
approximation introduced in the eikonal treatment of 
FSI, specifically the absence of any spin-orbit effect, 
does not affect the agreement between the solid lines 
and the data. Similarly, the sensitivity to the off-shell 
ambiguity in the electromagnetic current operator is 
relatively weak. After all, it is well known that the 
cross section is not particularly sensitive to the theo- 
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retical uncertainties in the description of FSI and the 
off-shell prescriptions within a range of about 10%. 
In Ref. |23|| more sensitive observables in the unpolar- 
ized cross section were considered. Good agreement 
with these data was maintained but, at the same time, 
the limitations of this approximation emerged, partic- 
ularly in the left-right asymmetry. Here, for the same 
kinematic conditions we extend the analysis to polar- 
ization observables. 
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FIG. 2: Polarization transfer components P' L ,P' T for the 
ie O(e, e'p) reaction at Q 2 = 0.8 (GeV/c) 2 in perpendicular 
kinematics M] leading to the 15 N p|, p| and s| residual 
states. Solid, dashed, dot-dashed, and dotted lines as in 
Fig. | 
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FIG. 3: The ratio P' T /P' L for the le O(e,e'p) reaction at 
Q 2 — 0.8 (GeV/c) 2 in perpendicular kinematics jL5| lead- 
ing to the 15 N pi, p| and s| residual states. Notations 
as in Fig. 



In Figs. |2| and |2| the polarization transfer compo- 
nents P' L ,P' T and their ratio P' T jP' L are shown as 
functions of the missing momentum p m , respectively, 
for the 16 0(e, e'p) reaction at Q 2 = 0.8 (GeV/c) 2 and 
constant (q, ui) for the transitions to the 15 N ground 
state p\, the first p| state at E m — 6.32 MeV and the 
weak peak with quantum numbers rising above a 
continuum background at E m ~ 28 MeV [}15| . Solid, 
dashed, dot-dashed, and dotted lines refer to the same 
calculations as in Fig. [j]. For these observables and at 
this kinematics, the sensitivity to off-shell effects is at 
most < 15%. The overall agreement with the data 
is still good, particularly for the microscopic calcula- 
tions with the p| quasihole state that performs even 
better than the RDWIA analysis presented here or ob- 
tained by other groups [jl5|. This fact is remarkable, 
since the RDWIA analysis depends on mean-field phe- 
nomenological potentials with several parameters fit- 
ted to the considered target and energy domain, while 
the calculation with quasihole states is basically pa- 
rameter free. In fact, from Eq. (Q) it is easy to verify 
that the polarization observables are given by ratios 
between a specific spin projection of the cross section 
and the unpolarized cross section, eliminating any sen- 
sitivity to the spectroscopic factor which is anyway 
fixed from the very beginning to the low-energy data 
of Ref. ]35|]. Moreover, the calculations of the solid 
lines include an attempt of a microscopic description 
of FSI in the framework of the eikonal approximation 
in a way which is consistent with the description of 
the bound state. The limitations of such an approach 
are more evident in the 3 = \ case, where, contrary 
to the RDWIA analysis, the absence of any spin-orbit 
effects is most likely responsible for the worse agree- 
ment. In any case, the second P' T data point for both 
pi and si shells appears not reproducible in both cal- 
culations, causing the theoretical ratio P' T jP' L to 
deviate substantially from the experiment. 

In summary, we have analyzed the unpolarized and 
polarized proton knockout reactions on 16 O at differ- 
ent kinematics with two theoretical approaches. The 
RDWIA is based on a relativistic mean-field descrip- 
tion of the proton bound state and on the effective 
Pauli reduction of the final Dirac spinor, leading to a 
Schrodinger-equivalent mean-field description of resid- 
ual FSI and to an effective relativistic electromagnetic 
current operator which depends on the Dirac scalar 
and vector potentials. The same kind of Pauli reduc- 
tion (and resulting current operator) is used for both 
initial and final states in the second approach, where 
a microscopic description of the bound state proper- 
ties is obtained by solving the Dyson equation with 
a nucleon self-energy which includes realistic short- 
range and tensor correlations for 16 O. As an attempt 
towards full consistency, the proton scattering wave is 
then generated in the eikonal approximation by micro- 
scopically calculating the damping of a plane wave as 
a solution of the Dyson equation for the nucleon self- 
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energy including the same realistic short-range and 
tensor correlations between the struck proton and the 
surrounding nucleons in nuclear matter. A systematic 
good agreement with data is observed for both unpo- 
larized and polarized reactions at low and high Q 2 
by using the same spectroscopic factors fixed by the 
low-Q 2 analysis, thus indicating that a high degree of 



internal consistency has been reached. 
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